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We calculate the density profiles and density correlation functions of the one-dimensional Bose gas 
in a harmonic trap, using the exact finite-temperature solutions for the uniform case, and applying 
a local density approximation. The results are valid for a trapping potential which is slowly varying 
relative to a correlation length. They allow a direct experimental test of the transition from the 
weak coupling Gross-Pitaevskii regime to the strong coupling, 'fermionic' Tonks-Girardeau regime. 
We also calculate the average two-particle correlation which characterizes the bulk properties of the 
sample, and find that it can be well approximated by the value of the local pair correlation in the 
trap center. 
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I. INTRODUCTION 

The simplest investigations into a many-body sys- 
tem like a Bose-Einstein condensate comprise studies 
of thermal equilibrium properties, and the physics of 
small fluctuations around thermal equilibrium. For one- 
dimensional systems, very similar behavior is found us- 
ing either photons in optical fibres or ultra-cold atoms in 
waveguides. Although techniques are not yet as exper- 
imentally advanced in the latter case, preliminary the- 
ory and some experimental measurements have already 
taken place. The atomic systems have the advantage 
that relatively long interaction times, large interaction 
strengths and low losses are possible, thus potentially 
allowing stringent tests of underlying quantum correla- 
tions. In this paper, we extend previous studies of cor- 
relations to include the experimentally realistic case of 
atoms in a waveguide with a harmonic longitudinal con- 
fining potential. The treatment is at finite temperature, 
and makes use of exact results for the uniform gas, to- 
gether with a local density approximation. 

For strong radial confinement, these types of system 
are examples of one-dimensional quantum (ID) gases 
IH El LJ They have the important property that in 
many cases their energy eigenstates are exactly solvable 
H,HBS EHBUllIli resulting in a greatly increased 
fundamental understanding of the relevant quantum field 
theory. For this reason, the study of ID systems plays 
an important role in the physics of quantum many-body 
systems. It is possible to make first-principle predictions 
without introducing added approximations like pertur- 
bation theory. This permits direct experimental tests of 
the underlying many-body quantum physics, as has been 
demonstrated in photonics with squeezed solitons in op- 
tical fibres fTi ]. 



For ultra-cold atomic systems with repulsive interac- 
tions, the most interesting and exciting feature is the 
predicted transition of an interacting gas of bosons to 
a 'fermionized' Tonks-Girardeau regime at large cou- 
pling strengths and low densities. 

We have recently made use of the known exact solu- 
tions to the uniform one-dimensional (ID) interacting 
Bose gas problem to calculate the exact local second- 
order correlation function at all densities and interaction 
strengths [TEl ITtj Il7j . This is the most direct indication 
of 'fermionic' behaviour, since this correlation function 
is strongly reduced at low density and strong coupling - 
similar to the case of fermions, where it vanishes exactly, 
due to the Pauli exclusion principle. 

The first experimental evidence of reduced or 'anti- 
bunched' correlations in a ID Bose gas has recently been 
demonstrated in Ref. [3| . However, current experiments 
typically take place in traps, with a longitudinal trap- 
ping potential. Provided the trap potential varies slowly, 
this environment is sufficiently close to a uniform one so 
that the exact solutions can still be used locally, in an 
approximation called the local density approximation. 

In this paper, we make use of the local density ap- 
proximation (LDA) to calculate the density profile and 
finite temperature local pair correlation function of a ID 
Bose gas trapped in a harmonic potential. The results 
are valid for sufficiently low longitudinal trap frequen- 
cies, and make use of the exact solutions to the plane- 
wave Lieb-Liniger model at finite temperature, to- 
gether with the Hellmann-Fcynman theorem. We mostly 
focus on regimes with quantum degeneracy. This requires 
temperatures T < Tq, where Tq = NHuj z is the temper- 
ature of quantum degeneracy of the trapped sample as a 
whole, N is the total number of particles, and uj z is the 
axial trap frequency. 
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Our main results show how 'fermionization' can be 
readily detected through a simple measurement of the 
pair correlation averaged over the trap. This is very close 
to the correlation at the trap center as predicted us- 
ing the Lieb-Liniger uniform model, and in principle can 
be measured via photoassociation of trapped atoms, or 
other related two-body inelastic processes whose rates 
are governed by the local pair correlations 01- I n a d- 
dition, an indirect measure of the pair correlations can 
be obtained [Til Il9j via the measurement of three-body 
recombination rates as recently demonstrated experimen- 
tally 0. 



II. ONE-DIMENSIONAL BOSE GAS 

One dimensional quantum field theories have the im- 
portant and useful property that they are often exactly 
solvable. This is not generally the case for higher- 
dimensional quantum field theories. Thus, the study of 
these ID models can lead to an insight into the nature of 
quantum field theory for interacting particles, that is not 
possible from the usual perturbative approaches. In this 
section, we review the physics of these exact solutions for 
interacting bosons in the uniform ID case, and introduce 
the theoretical framework for treating a non-uniform gas 
within the local density approximation. 



A. Hamiltonian 

The study of exact solutions for the one-dimensional 
Bose gas started with Girardeau's seminal work on 
hard-core, or impenetrable bosons. In this model, there 
is a remarkable and exact correspondence between the 
measurable correlation functions of free fermions, and 
those of strongly interacting bosons. In the ID Bose 
gas model with a delta-function interaction, solved by 
Lieb and Liniger , the particles can pass through each 
other, so they are no longer impenetrable. This provides 
a realistic description of a wave-guide with transverse di- 
mensions larger than the 'core' of a particle in the wave- 
guide. Under these circumstances, there may only be 
a single relevant transverse mode, yet particles are able 
to exchange their positions as they propagate past each 
other. 

Thus, we start by reviewing the theory of a gas of N 
bosons interacting via a delta- function potential in one 
dimension. The ID Bose gas has a short-range repulsive 
interaction between particles which is characterized by 
just one coupling constant. In second quantization, the 
Hamiltonian is 

H = — [ dzd z &d z V+~ f dz&&^ 
2m J 2 J 

dzV(z)&&, (2.1) 



where ^(z) is the bosonic field operator, m is the atom 
mass, g > is the coupling constant, and V(z) is the 
trapping potential which we assume is harmonic with 
V(z) — muj 2 z 2 /2, while uj z is the trap oscillation fre- 
quency in the axial direction. To treat the uniform gas 
we set V(z) = 0. 

For Bose gases in highly elongated cylindrical traps 
(uj z <C lu± , where tu± is the frequency of the transverse 
harmonic potential) such that the sample can be de- 
scribed by the above ID model, the coupling constant g 
is expressed through the 3D scattering length a |2fj( . For 
a positive scattering length a which is much smaller than 
the amplitude of transverse [x, y-direction) zero point os- 
cillations, or the transverse harmonic oscillator length, 



one has 



2h 2 a 



fj ■ 



muj±, 



2hhj i a. 



(2.2) 



(2.3) 



The ID regime is reached if lj_ is much smaller than the 
thermal de Broglie wavelength At = (2ttH 2 /mT) 1 / 2 and 
a characteristic length scale l c |21| responsible for short- 
range correlations. On the same grounds as at T = [lfij . 
one finds that for fulfilling this requirement it is sufficient 
to satisfy the inequalities 



a < l ± < {l/n(0), A T }, 



(2.4) 



where u(0) = ^\E ft (0)"I f (0)^ is the ID (linear) density in 
the center of the trap, z = 0. 

B. Ground-state solution for the uniform gas 

We now give a brief overview of the uniform (V(z) = 0) 
Bose gas problem describing a gas of N bosons interact- 
ing via a pair-wise repulsive delta-function potential in 
a ID box of length L with periodic boundary condition 
0- In the thermodynamic limit (N, L — > oo, while the 
ID linear density n = N/L is kept constant), the solu- 
tion to the energy eigenstates is found |7j using the Bethe 
ansatz 22]. In this solution, all relative wave- functions 
are assumed to have a plane- wave form - except for finite 
changes in gradient at each collision where two particle 
coordinates are equal: 



iV 



*">= K^ri' 



n i 



3>i 



-sgn(zi - Zj) 



x ¥(z 1 )¥(z 2 )...¥(z N )\0) 



(2.5) 



Here, we have used units in which h = 2m = 1, while 
introducing the Lieb-Liniger notation Q of c = mg/h 2 . 
Also, sgn(zi — Zj) is the sign function, and fcj is the 
"quasi-momentum" . The quasi-momenta are determined 
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from the delta-function slope-change requirements, that 
is, fcj+i — kj is determined by the boundary conditions 
at Zi — Zj. In the limit of a large sample, and defining 
fej+i — kj = l/[L/(fcj)], one can approximate f(kj) by a 
continuous function f(k) which is the density of quasi- 
momenta. The distribution of quasi-momenta is then ob- 
tained as the solution to the following integral equation: 



so that 



27r/(fc) = 1+ / K(k-p)f(p)dp. 
Here, the kernel function K{k) is given by 
K{k)= -'' 



k 2 ' 



(2.6) 



(2.7) 



and kp is the maximum quasi- momentum which deter- 
mines the particle number density n = N / L via 



f(k)dk 



(2.8) 



— kp 



The corresponding ground-state energy is given by 

rkp 

'2 



E = L 



f(k)k 2 dk, 



(2.9) 



and is often written as Eq = Nn 2 e(j), being an implicit 
function of n, via a dimensionless function e(j) of the 
parameter 7 = c/n. 

Restoring the physical units, this gives an energy per 
particle: 



Eo/N = ^n 2 e( 7 ), 



where 



7 



mg 
h 2 n 



(2.10) 



(2.11) 



The dimensionless parameter 7 which characterizes the 
strength of interactions is in fact the only parameter 
needed to describe the uniform ID Bose gas at zero 
temperature. The limit of 7 < 1 corresponds to the 
weakly interacting Gross-Pitaevskii (GP) regime, where 
the mean-field Bogoliubov theory works well. The oppo- 
site limit of 7 3> 1 corresponds to the strongly interacting 
or Tonks-Girardeau (TG) regime, and as 7 — * 00 one re- 
gains Girardeau's results for impenetrable bosons. 

The solution to the ground state energy Eq can be 
used, together with the Hellmann-Feynman theorem |23| , 
for calculating an important observable - the normalized 
local pair correlation 



(2.12) 



The pair correlation is found by taking the derivative 
of the ground state energy with respect to the coupling 
constant g, owing to the fact that 



2 (2) (0) 



de(r/) 



(2.14) 



The pair correlation </ 2 )(0) for the zero temperature 
uniform ID Bose gas has been calculated using the Lieb- 
Liniger exact solution for e( 7 ). The results are given 
in Ref. [lj|. Here, we will extend these results (see Sec- 
tion III) to the case of a trapped (non-uniform) Bose 
gas using the local density approximation, and to finite 
temperatures as well. 



C. Uniform gas at finite temperature 

The excited states of the uniform ID Bose gas can be 
calculated in a similar way, with each excited state cor- 
responding to the removal of a quasi-momentum with 
|fc| < kp - called a hole - and the creation of a quasi- 
momentum with \k\ > kp. In 1969, C.N. Yang and CP. 
Yang @ worked out the finite-temperature density ma- 
trix solution for the Lieb-Liniger model, by constructing 
the free energy and taking into account the entropy of 
all the different excited states. This was used in a sub- 
sequent work to calculate numerically the pressure of 
the function of temperature. 

At thermal equilibrium, we now assume that the den- 
sity of quasi-momenta f(k) has no upper cut-off, and 
that it consists of two types of terms - occupied quasi- 
momenta [with density f p (k)) and unoccupied or 'hole' 
quasi-momenta [with density fh{k)], i.e. f(k) — f p (k) + 
fh(k). The overall integral equation now has the form 
(using units in which h = 2m = 1 and c = rag/h 2 ) 

2irf(k) = 2n [f p (k) + f h (k)} 

= 1+/ K(k-p)f p (p)dp. (2.15) 



Hence, the particle density n is obtained from the oc- 
cupied or particle quasi-momenta 

f P (k)dk, (2.16) 

while the total energy is now 

/oc 
f p (k)k 2 dk. (2.17) 
-00 

However, there is also an entropy involved, since there 
are many wavefunctions that are nearly the same, within 
a given range of values of f p (k) and fh{k). In fact, the 
number of choices compatible with a given dk value is 



[f(k)Ldk]\ 



[f p (k)Ldk]\[f h (k)Ldk]\ 
Thus, the entropy is 



(2.18) 




L 



(^(z)^( Z )*(z)*(z)), (2.13) S = L T [f\nf-f p \nf p -f h \nf h ]dk. (2.19) 

J — OO 
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Minimizing the total free energy F = E — TS gives the 
thermal equilibrium distribution of holes and particles, 
where we choose temperature to be in energy units, so 
that fcs = 1. The minimization at a fixed average particle 
number requires the use of a Lagrange multiplier /i, and 
gives the result that the distribution fp(k) satisfies the 
integral equation 



27T/ p (fc) 



1 



a e{k)/T 



= 1- 



K(k-p)f p (p)dp, (2.20) 



where the excitation spectrum e{k) is calculated from a 
second integral equation 



e{k) 



— J K(k-p)ki(l + e- e W/ T )dp. (2.21) 



Here, \x can be shown to coincide with the chemical po- 
tential of the system, while the entropy and the free en- 
ergy per particle are found from: 



S/N = - 
n 



1 



/(fc)ln(l + e- £(fc)/T 
f p {k)e(k)dk, 



dk 



(2.22) 



FIN = u — — (°° In fl + e- £(k)/T ) dk. (2.23) 

In addition, using the thermodynamic identity F = 
—PL + /iN, one can arrive at the following simple result 
for the pressure of the gas 

P(u, T) — / In ( 1 + e - £ ( fe )/ T ) dk. (2.24) 

To calculate the pair correlation g( 2 ' (0) for a finite tem- 
perature gas one can again use the Hellmann-Feynman 
theorem [23| . Here, we consider the canonical partition 
function Z = exp(-F/T) = Trexp(--ff/T), where the 
trace is over the states of the system with a fixed particle 
number N, at temperature T. Taking the derivative of 
F = — Tlog Z with respect to the coupling constant g we 
obtain 



dF 1 



-H/T 



dH 
dg 



(2.25) 

Introducing the free energy per partic le f — F/N and 
restoring the physical units, this gives 



;lef 

m 



3 (2) (0)- 



Ln 2 



dF 
dg 



h 2 n V 



N .T " " \ 9 1 

(2.26) 

Here, r = T/Td is a dimensionless temperature parame- 
ter, with Td = h 2 n 2 /2m being the temperature of quan- 
tum degeneracy for a uniform gas. Hence, we have 



The pair of the dimensionless parameters 7 and r com- 
pletely characterize the properties of a finite temperature 
uniform gas. 

Alternatively, the local pair correlation g^> (0) can be 
calculated within the grand canonical formalism. Here, 
we consider the grand canonical partition function Z = 
exp(-0/T) = Trexp[(iiN-H)/T}, where Q = F—fiN = 
—PL is the grand canonical thermodynamic potential 
and P is the pressure. The trace is over the states of 
the system, at a fixed chemical potential jj, and temper- 
ature T . Taking the derivative of £1 = —Tlog Z with 
respect to the coupling g we obtain: 

dil d log Z 

dg dg 

= 4lv{«^)e xp[ (^-F)/T]} 



- (*t(z)*t(*0*(z)*(z)) 



(2.28) 



Thus, the normalized pair correlation g^(0) can be 
calculated using 



/ 2) (0) = 



Ln 2 



air 

9 9/ p 1 



1 



dP" 
d 9 ; n 1 



(2.29) 



This requires the use of Eq. (|2.24|) for the pressure, 
which in turn is found after solving the Yang- Yang in- 
tegral equations I2.20|l and (|2.21|) . 

The local pair correlation for a finite temperature uni- 
form gas has been first calculated in Ref. [l|j using 
the exact solutions to the Yang- Yang integral equations 
(|2~2"0"|) and i(2~2"T|l . together with Eqs. l(2~2lj) and l(2~27)j> . 
In Sections IV- VII, we will extend these results to the 
case of a trapped gas using the local density approxima- 
tion. 



D. Quasi- uniform approximation 

In a quasi-uniform approximation, we suppose that 
the system can be divided into small regions of size Az 
which is larger than a characteristic short-range correla- 
tion length l c . In each of this regions we assume that the 
inhomogeneity of the gas is negligible so that it can be 
treated as a uniform gas. 

In this case, the trapping potential V(z) is replaced by 
a step-like function V(zj) = mu) 2 z 2 /2 which is constant 
within each region from Zj to and undergoes step- 
like changes at the boundaries of the adjacent regions. 
Here, the size Az takes the role of the length L from the 
Yang- Yang solution that applies to each region. 

We now consider an ansatz in which the overall density 
matrix has the structure of an outer-product of canonical 
solutions, with Nj being the average number of particles 
in the j-th region: 



n 2 . 



>N nt 



(2.30) 
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Next we look for an approximate solution in which the 
effective Hamiltonian is assumed to introduce no cou- 
pling between the regions. To obtain this we must now 
minimize the total free energy given by 



(2.31) 



3=1 



This requires us to include a constraint on the total par- 
ticle number: 



N 



E^ 

3=1 



(2.32) 



Hence, it is appropriate to use a Lagrangian formulation 
with 



c = E & - TS > - i* N i 



(2.33) 



3 = 1 



We note here that the Lagrangian C is now simply 
a sum over independent regions, with each term corre- 
sponding to that for a single uniform Bose gas. As we 
are only constraining the total particle number, not the 
number in each region, the chemical potential is the same 
for each term. Since there is no explicit coupling be- 
tween the regions, the Lagrangian is minimized when we 
satisfy the Yang- Yang equations in each separate region, 
but with the same (global) chemical potential fia at all 
locations. 



E. Local density approximation (LDA) 

In more detail, we have shown that for a large system, 
where the density profile varies in a smooth way, the 
system behaves locally as a piece of a uniform gas. This 
can be described locally as a uniform gas with chemical 
potential equal to the local effective chemical potential 



fi(z) = Mo - V(z) = no 



1 



2 2 



(2.34) 



where /^o is the global equilibrium chemical potential. 

For the LDA to be valid, the short-range correlation 
length l c (z) should be much smaller than the character- 
istic inhomogeneity length k n h{z)- These length scales 
depend on the displacement from the trap center z, and 
the LDA validity criterion reads: 



l c (z) < hnh(z) 



n{z) 



\dn(z) / dz\ 



(2.35) 



The short-range correlation length l c (z) is defined lo- 
cally via the density distribution n(z). At low temper- 
atures, the correlation length l c {z) can in general (irre- 
spective of the interaction strength) be expressed via the 
local chemical potential fj,(z): 



lc(z) 



y/mn{z)' 



(T « T Q ). 



(2.36) 



In the weakly interacting Gross-Pitaevskii (GP) regime 
the relation between the chemical potential /i(z) and 
the density n(z) is n(z) = gn(z), and we obtain that 
the correlation length coincides with the healing length 
l c (z) = h/iJmgn(z). In the strongly interacting Tonks- 
Girardeau (TG) regime one has (i(z) = TT 2 h 2 n 2 (z)/(2m), 
so that l c (z) ~ l/n(z), neglecting the numerical factor 
of order 1. At high temperatures, l c {z) is of the order of 
the thermal de Broglie wavelength At- 

The condition (|2.35l) is sufficient for using the LDA 
for calculating the density profiles and local correlation 
functions. The reason is that these correlations, in partic- 
ular the two-particle correlation g( 2 >(0), are determined 
by the contribution of excitations which have energies 
of the order of the chemical potential and wavelength of 
the order of l c . However, this is not the case for all cor- 
relation functions. For example, calculation of the finite- 
temperature single-particle correlation function would re- 
quire a strong LDA condition in which the sample size 
was much larger than the phase correlation distance [27j ■ 
In this sense, one may call Eq. (|2.35|l the 'weak' LDA 
criterion. However, within the LDA, no correlation func- 
tion can be calculated reliably over distance scales that 
are comparable to the sample size. 

Thus, the 'weak' criterion of validity of the LDA re- 
quires that variations of the density occur on a length 
scale that is much larger than l c {z), in which case the 
gas is treated locally as a piece of a uniform gas. From 
the definition of l% n h{z), one can easily see that the LDA 
is easier to satisfy in the center of the trap where the 
density profile is almost flat than near the tails of the 
distribution where the density drops rapidly. However, 
for measurements that average over an entire trap, it is 
the central region that plays the most important role. 



III. ZERO-TEMPERATURE TRAPPED GAS 
A. LDA criterion at T = 

Here, we analyze the implications of the LDA criterion 
(|2.35l) for a zero temperature gas. At T = 0, a uniform 
ID Bose gas can be characterized by a single dimension- 
less interaction parameter 7, Eq. ( I2.11|l . Depending 
on its value, one has two well-known and physically dis- 
tinct regimes of quantum degeneracy. For 7 <C 1, i.e. 
at weak couplings or high densities, the gas is in a co- 
herent or Gross-Pitaevskii (GP) regime. In this regime, 
long-range order is destroyed by long-wavelength phase 
fluctuations 0> HE HE 113 an d the equilibrium state 
is a quasi-condensate characterized by suppressed den- 
sity fluctuations. For strong couplings or low densities, 
7> 1, the gas reaches the strongly interacting or Tonks- 
Girardeau (TG) regime and undergoes "fermionization" 
[E 0|- The term "fermionization" is used here in the 
sense that the wave function strongly decreases as parti- 
cles approach each other. 

For a trapped (non-uniform) gas one can introduce a 
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local interaction parameter 

70) 



mg 
h 2 n(z)' 



(3.1) 



which changes with the density distribution n(z) and can 
be used for characterizing the local properties of the gas. 

From the definition of j(z) it is clear that as one moves 
from the center of the trap towards the tails of the density 
distribution where n(z) — > 0, the gas either enters the 
TG regime where -f(z) ^>> 1, or else the LDA itself breaks 
down. 

Moreover, in the T = case the Lieb-Liniger solu- 
tion within the LDA gives a density profile that van- 
ishes \ n(z) = 0] beyond a certain distance R from the 
origin [28l |2|j. This distance is called the Thomas- 
Fermi radius and it is determined from the condition 
/i(-R) = Mo — mu 2 R 2 /2 = which gives: 



R 



( 2^o 
\mu> 2 



1/2 



(3.2) 



Since n(z) vanishes exactly at \z\ > R |30j , it is clear 
that the LDA criterion (|2.35l) can only be satisfied up to 
a certain maximum distance from the trap center \z\ ~ 
R — Sz, displaced from R by 5z (Sz -C R). We would 
like therefore to determine the displacement Sz such that 
the LDA is valid for < \z\ < R — 5z and breaks down 
beyond \z\ ^ R — Sz. 

As we are interested in calculating the density pro- 
files and the local two-particle correlation function, wc 
will only focus on the 'weak' LDA criterion, Eq. H2.35fl . 
First, we rewrite the inhomogeneity length scale from Eq. 
13 in the following equivalent form 



knh(z) = n(z) 



dfi(z) / dn(z) 



dfj,(z)/dz 



(3.3) 



Using the explicit expression /i(z) ■ 
taking its derivative, we obtain at \z\ 



fiQ — mLo^z /2 and 
- R-Sz (Sz < R): 



linh (z*) 



n(z) 
muj 2 R 



dfx(z) 



dn(z) 



(3.4) 



Combining Eqs. (|2.36|) and 13.4|l . one can rewrite the 
LDA criterion l c (z) <SC hnh(z) in the following equivalent 
form (again neglecting numerical factors of order one): 



/i(z)\ 3 ^ 2 no d\nji(z) 



Mo 



ftw d\un(z) 



> 1. 



(3.5) 



Next, we note that in the limiting GP and TG regimes 
the derivative dlnfj,(z)/dlnn(z) is equal, respectively, to 
one and two, so in general 1 < dhifj,(z)/dhin(z) < 2. 
Therefore, the exact numerical value of this quantity can 
be replaced by unity in all regimes, and the LDA criterion 
becomes 



mO) 

Mo 



3/2 



MO 



» 1. 



(3.6) 



Finally, expanding fi(z) near the edge of the cloud 
where \z\ ~ R — Sz (Sz <iJ), we obtain that fi(z)/[iQ ~ 
Sz/R, so that the criterion of applicability of the LDA is 
reduced to a simple requirement: 



Sz 
~R 



> 



hjj z 
Mo 



2/3 



(3.7) 



For 7(0) <C 1, the local value of j(z) at z ~ R — Sz is 



7 (z) 



mgR 
H 2 n(0)Sz 



We sec that if 7(0) < (fiw z /> ) 2/3 then j(z) < 1 at 
z ~ R — Sz, implying that the gas stays in the GP regime 
at all locations z until the LDA breaks down. If, on 
the other hand, 7(0) ^> (fvjj z / Mo) 2 ^ 3 then 7(2) ^> 1 at 
z ~ R—Sz so that the gas first approaches the TG regime 
and then the LDA breaks down. 

For 7(0) 3> 1 the gas is in the TG regime at all loca- 
tions z, until the LDA breaks down. This is because "f(z) 
is always larger than 7(0), and hence 7(2) ^> 1. 

In the limiting cases of 7(0) <C 1 and 7(0) 3> 1, 
the LDA criterion l|3.7|l can be conveniently rewritten 
in terms of 7(0) and the total number of particles TV. 
In doing so we use the fact that the chemical poten- 
tial is given by fio — gn(0) for 7(0) <C 1, and by 
^0 = Tr 2 h 2 n 2 (0)/(2m) for 7(0) 3> 1 ■ In addition, we use 
the known relationship between the peak density n(0) 
and N in each case [see Eqs. (I3.14f) and (13.181 ) below]. 
As a result, we obtain that the LDA criterion (|3.7() can 
be rewritten as follows, in the GP and TG regimes, re- 
spectively: 



2/3 



(3.8) 



Sz 1 fku z \ 

>> vRoyJ 



2/3 



R 

Sz 

"z7 



> 



2mui z 



ir 2 H 2 n 2 (0) 



~ 7 ( )i/3/v2/3 

2/3 1 



N 2 / 3 '' 



, [ 7 (0) «1], (3.9) 



[ 7 (0) »1]. (3.10) 



As we see, for any small but finite 7(0) in the GP 
regime, the right hand side of Eq. 13.9fl can be made 
small by increasing the total number of particles N. For 
a fixed coupling g and a constant density n(0) [such that 
7(0) stays constant], the increase of the particle num- 
ber N has to be accompanied by a reduction of the trap 
frequency oj z . Thus, the ratio Sz/R can also be made 
small, so that that the LDA criterion in the GP regime 
[7(0) <C 1] is satisfied for almost the entire sample, up 
to the location z = R — Sz very close to the edge of the 
cloud. Similar considerations apply to the TG regime 
[7(0) 3> 1], where the requirement on the (large) total 
number of particles N is less stringent than in the GP 
regime. 



7 



B. Pair correlations at T = 



Here we discuss the local pair correlation 



g^(z,z)^ 



# t (z)* t (z)*(z)*(z) 
n 2 (z) 



(3.11) 



in a zero temperature trapped gas within the LDA. The 
calculations are done using the solution to the Lieb- 
Liniger equation 1|2.6|) and the Hellmann-Feynman the- 
orem HH , Eq. ifTTIjl . Here, g {2) (0) is now replaced by 
g^(z, z) and 7 is to be understood as the local value of 
j(z). Thus, to calculate g^ 2 \z,z) as a function of the 
distance from the trap centre, one can use the uniform 
results in which the interaction parameter 7(2) is found 
from the density profile n(z), for different values of 7(0). 

The implementation of the LDA, using the local effec- 
tive chemical potential n{z), Eq. (|2.34|) . is carried out 
by means of first calculating the chemical potential (i as 
a function of n, and then inverting this dependence for 
obtaining n(z). This gives n(z) [and hence 7(2)] as a 
function of for a given value of the interaction pa- 
rameter 7(0) at the trap center. 

Depending on the value of the coordinate-dependent 
interaction parameter 7(2), Eq. I|3.1|l . we have the fol- 
lowing limiting behavior of the pair correlation function. 

In the Gross-Pitaevskii (GP) limit of a weakly inter- 
acting gas, the pair correlation in the uniform case is 
g(2) ~ 1 _ 2^/7/77, 7 <C 1 0]. For a trapped gas, replac- 
ing 7 by 7(2) gives: 



g {2) {z,z)^i 



7(0) 



1 -z 2 /R 2 ' 



7 (z)<l. (3.12) 



where we have used the relationship between 7(2;) and 
n(z) and the fact that in the GP regime the density pro- 
file is given by the familiar Thomas-Fermi parabola: 

n(z) = n(0)(l - z 2 /R 2 ), (3.13) 

and n{z) = for \z\ > R. Here, the peak density n(0) 
and the radius R are given by: 



n(0) 



R 



V 32 5 
/ 3Ng \ 1/3 

\ llXlLOt 



(3.14) 



(3.15) 



In the Tonks-Girardeau (TG) limit of strong inter- 
actions, the uniform gas pair correlation is g^ ~ 
47r 2 /(37 2 ), 7 >1 [15|. In the trapped gas case, replacing 
7 by 7(2) gives 



9 {2) {z,z) 



4tt 2 (1- z 2 /R 2 ) 

WW) : 



7 (z) » 1, (3.16) 



where we again used the relationship between 7(2) and 
n(z) and the fact that the density profile in the TG 
regime is given by the square root of parabola: 




FIG. 1: Examples of the density profiles n(z)/n(0) of a zero 
temperature ID Bose gas in a harmonic trap as a function 
of the dimensionless coordinate z/R, for different values of 
the interaction parameter 7(0). The solid lines are the results 
of the exact numerical solution of the Lieb-Liniger equations 
within the LDA. The dashed and the dotted lines are the 
analytic results given by the Thomas-Fermi parabola in the 
GP regime and the square root of parabola in the TG regime, 
respectively. 



and n(z) = for \z\ > R. Here, the peak density n(0) 
and the radius R are 



n(0) 
R 



f2mNuj z \ 1/2 



\ ir 2 h 
2hN\ 1/2 



muj 



■ 



(3.18) 
(3.19) 



n(z) = n(0)(l-z 2 /i? 2 ) 1/2 , 



(3.17) 



In Fig. ^ we present the density profiles n(z) as a 
function of the dimensionless coordinate z/R, for differ- 
ent values of 7(0). The full lines represent the results of 
the numerical calculation within the LDA, which repro- 
duce the results of Ref. [2^], while the dashed and the 
dotted lines represent, respectively, the above analytic 
results in the GP and TG regimes. 

Figure shows the local pair correlation g^(z,z) as 
a function of z/R, for different values of 7(0). As we 
see, for the case of weak interactions 7(0) <§; 1, the pair 
correlation is close to unity in the central bulk part of 
the distribution. This is an expected result for the co- 
herent or GP regime. As one approaches the tails of the 
distribution, where the density is reduced and 7(2) be- 
comes larger than one, the gas locally approaches the TG 
regime of 'fermionization'. Here, the pair correlation is 
suppressed below the coherent level, g^ 2 \z,z) < 1. 

For the cases where 7(0) > 1, including the TG regime 
of 7(0) S> 1, the pair correlation is suppressed across the 
entire sample. In the limit 7(0) — » 00, the gas acquires 
pronounced "fermionic" properties so that the wave func- 
tion strongly decreases as particles approach each other, 
thus resulting in g^(z, z) — > 0. 
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X(0)=0.01 




0.4 0.6 

z/R 



FIG. 2: The local pair correlation of a trapped ID Bose gas at 
zero temperature, </ 2 '(z, z), as a function of the displacement 
from the trap center z/R, for different values of 7(0). The 
full lines are the results of numerical calculation, while the 
dashed lines for 7(0) = 0.01,0.1 and 7 (0) = 10, 1 00 are the 
respective analytic results of Eqs. 13.12H and 13. 1611 shown for 
comparison. 



IV. FINITE TEMPERATURE TRAPPED GAS 
A. Key parameters 

An obvious choice of dimensionless interaction and 
temperature parameters for describing a trapped ID 
Bose gas within the LDA consists in using the local value 
of the interaction parameter 7(0) and the reduced tem- 
perature t(0) in the trap center. 

These are the same parameters that are used in the uni- 
form gas treatment [0|, Eqs. i|2.11|) and l|2.27[l . except 
that now we define them via the local density n(z). Thus, 
in general, we define the local interaction parameter 7(2:) 
and the local reduced temperature t(z), according to: 



7O) 
r(z) 



h 2 n(z) ' 
T 2mT 



T d (z) h 2 n 2 (zy 



(4.1) 
(4.2) 



conditions where the mean inter-particle separation be- 
comes of the order of the thermal de Broglie wavelength. 

The values of these parameters at the trap center, 7(0) 
and t(0), completely characterize all relevant properties 
of the gas within the LDA, including the associated den- 
sity profiles n(z), the resulting total number of particles 
N, as well as the correlation functions and the thermo- 
dynamic properties. 

A completely equivalent pair of the interaction and 
temperature parameters, which is, however, more suit- 
able for practical purposes is the local value of 7(2) and 
a new temperature parameter t defined via: 



r(0) _ t{z) 



T 



7 2 (0) 7 2 ( 2 ) mg 2 /(2h 2 )' 



(4.3) 



According to this definition, the temperature is measured 
in units of the characteristic energy E{, = mg 2 /(2h 2 ). 

The advantage of using t as the dimensionless temper- 
ature parameter is that it is independent of the density 
and gives a direct measure of the global temperature of 
the gas, which in equilibrium is the same for the entire 
sample. This allows us to easily explore the "interaction- 
temperature" parameter space [7(0) — t] in a systematic 
way. For example, considering different values of 7(0) 
while t is kept constant would correspond to physical 
conditions under which the peak density of the gas ro(0) 
is varied while the absolute temperature T is kept un- 
changed. A novel experimental technique that imple- 
mented this approach for achieving a Bose-Einstein con- 
densation in a 3-D gas has recently been demonstrated in 
Ref. H3. 

Other alternative choices are possible for characteriz- 
ing the interactions and temperature of a trapped gas in 
dimensionless units. For example, to characterize the sys- 
tem at different temperatures T while the total number of 
particles N is kept constant, one can define an alternative 
pair of global parameters which are more suitable for this 
case (see Sec. VI). Here, the global temperature parame- 
ter can be defined as — T/Tq, where Tq = Nhuj z is the 
global temperature of quantum degeneracy of a trapped 
gas (in energy units, ks — 1, where fcs is the Boltzmann 
constant). Irrespective of the interaction strength, a har- 
monically trapped Bose gas at T ^> Tq obeys the clas- 
sical Boltzmann statistics, whereas for T < Tq quantum 
statistical effects become important. This is clearly seen 
in the limit of a trapped ideal gas (7(0) — > 0) and in the 
opposite limit of a strongly interacting gas (7(0) — > 00). 
In the latter case, the problem maps onto the trapped 
gas of non-interacting fermions |(| . So, in both limits Tq 
appears explicitly as the temperature of quantum degen- 
eracy for the trapped sample as a whole. 



B. LDA criterion at finite T 



where T d (z) — h 2 n 2 (z)/(2m) is the local temperature 
of quantum degeneracy that corresponds (locally) to the 



Here we analyze the local density approximation for a 
finite temperature gas, and obtain simple criteria for its 
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validity in the limiting cases of very high and very low 
temperatures, T ^> Tq and T <C Tq. 

In the high temperature limit, T 3> Tq [in which case 
t(0) ~ Att(T/Tq) 2 ], the density distribution n(z) can be 
approximated by a Gaussian profile in all regimes, as the 
interaction between the particles is negligible compared 
to their thermal kinetic energies. For N particles at tem- 
perature T in a harmonic trap of frequency uj z , the den- 
sity profile is determined by the thermal distribution for 
a classical ideal gas described by Boltzmann statistics: 

n T (z) = -JL- exp(-z 2 /i4), (4.4) 

where the radius Rt characterizes the width of the Gaus- 
sian and is given by 

/ 2T 

R T = \ 2- ( 4 - 5 ) 

y mu) z 

In this high-temperature limit, the correlation length 
is given by the thermal de Broglie wavelength Ay , so that 
the LDA criterion (|2.35|1 gives: 

T ( 2T \ 1 ^ 2 T 

2 = tT7 R *> ( T » T Q)- ( 4 - 6 ) 

Since T Tq — Nftio z implies that T 3> Hcu z , the above 
LDA criterion can be satisfied for all locations z from the 
trap center up to distances equal to several characteristic 
widths Rt- For sufficiently large total number of parti- 
cles N, the ratio T/fuv z will be even larger so that the 
LDA will be valid for even larger distance from the trap 
center. 

In the opposite limit of low temperatures T -C Tq 
[in which case r(0) ~ 8a/2 7 (0)T/(3T q ) for 7 (0) < 1, 
and r(0) ~ it 2 T/T q for 7 (0) > 1], the density profile 
n(z) can be approximated by two contributions. The 
first one is for the central bulk part which will be close 
to the T = density profile up to a certain distance 
\z\ ~ R — 6z (Sz <C R) from the trap center. Here R 
is the zero-temperature Thomas- Fermi radius, Eq. I|3.2|l . 
The second contribution is for the tails of the distribution 
which can be approximated by a thermal Gaussian. 

As before, we will focus on the 'weak' LDA condition 
(|2.35|) involving the correlation length l c (z). For the cen- 
tral part of the sample, up to distances \z\ ~ R — Sz 
(Sz <C R), we can use the LDA criterion derived for 
the zero temperature gas, Eq. I|3.7|) . In the weakly and 
strongly interacting limits this can be rewritten - as be- 
fore - in terms of the interaction parameter 7 (0) and the 
total number of particles N [see Eqs. (13 .9|) and (|3.10[l ]. 

For the Gaussian tails of the distribution, i.e. at dis- 
tances \z\ > R, where the local correlation length is given 
by At, we use the above high-temperature result, Eq. 
(|4.6|) , and rewrite it in terms of the Thomas- Fermi radius 
R and the global chemical potential /iq. As a result, the 
LDA criterion for the Gaussian tails of a low-temperature 



gas can be written in the following form: 

/ J, \ 1/2 / V 1/2 

In the GP regime [ 7 (0) <C 1], where fio — 9 n (0), this 
gives: 

z/R<£t(0) 1/2 N. (4.8) 

Thus, in order that the LDA works in the tails of the den- 
sity distribution (\z\ > R) in the GP regime, one has to 
have t(0) 1 / 2 A^ ^> 1. This requirement can always be sat- 
isfied with a sufficiently large number of particles N. For 
example, for t(0) = 3.8 x 10 -3 [which can be obtained, 
for example, with 7 (0) = 0.01 and T/T Q = 0.01, accord- 
ing to the relationship r(0) ~ 8^2 7 (0)77(3T Q ) valid in 
this regime] one would need to have N > 870 in order to 
satisfy the LDA criterion for distances R < z <C 20R. 

In the TG regime [ 7 (0) ^> 1], where fio ~ 
Tr 2 h 2 n 2 (0)/(2m), Eq. I|4.7|l again reduces to the condi- 
tion given by Eq. (|4.8|) . Thus, in the TG regime the 
validity of the LDA in the tails of the distribution again 
requires that r(0) 1 / 2 A r ^> 1. However, now we have 
t(0)~t: 2 T/Tq, for T«T q . 

To summarize, in the low temperature limit the LDA 
criterion can be easily satisfied for the central bulk part 
of the density distribution and for the Gaussian tails. 
This leaves the question of validity of the LDA in the 
low-density region near z = R, where the density may 
vary more rapidly. 

We note, however, that at small finite temperatures 
the variation of the density profile around z = R is more 
smooth than in the T = limit, so that the LDA crite- 
rion may still be satisfied in this region, in contrast to the 
T = case where the LDA necessarily breaks down as one 
approaches the edge of the cloud at z = R. More impor- 
tantly, the LDA becomes valid again for distances past 
the small critical region around z = R, i.e. in the tails 
of the density distribution. This means that the results 
of calculation of the pair correlation function g^(z,z) 
at small finite temperatures should be valid everywhere 
except in a small region around z = R. At high tempera- 
tures the LDA criterion becomes less restrictive, and can 
be satisfied for the entire sample. 



C. Calculating the local pair correlation and 
density profiles 

The local pair correlation g( 2 \z,z), Eq. I|3.11(l . as a 
function of the location z from the trap center is calcu- 
lated using Eq. (|2.29H in which /i is replaced by the local 
chemical potential fi(z) — /iq — V(z) and where n is the 
local density n(z). The calculation is based on iterative 
numerical solution of the Yang- Yang exact integral equa- 
tions |8( for the excitation spectrum and for the distribu- 
tion function of "quasi-momenta" , l|2.20|l and (|2.21|) . For 
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a given set of values of (i(z), T and g, this gives the re- 
sulting density profile n(z), Eq. H2.16fl . and the pressure 
P, Eq. H2.24JI . Differentiating P with respect to g gives 
the local pair correlation g^ 2 \z,z). 

A convenient way to implement the numerical algo- 
rithm for solving the Yang- Yang equations is via a di- 
mensionless coordinate 



z 

Rt 



(4.9) 



where the length scale Rt is the thermal width of the 
classical Gaussian distribution t%t{z), given by Eq. I|4.5J) . 

Using the dimensionless coordinate £, the local chem- 
ical potential can be rewritten as 



h{RtO = mo - re 



(4.10) 



After setting up a lattice of £ values, {&}, the solution 
to the Yang- Yang equations proceeds as in the case of a 
uniform gas, with the input parameters being an array of 
the values of the local chemical potential ^ = /j,(Rx£i), 
the temperature T, and the coupling parameter <?, as 
described above. 

The final numerical results are then presented in terms 
of the dimensionless parameters 7(0) and the tempera- 
ture parameter r(0) (or i), where we note that 



n(z) = 7 (0) 
n(0) i(z)- 



(4.11) 



This makes the output results scalable with respect to 
the physical parameters, rather than dependent of their 
absolute values. 

The total number of particles in the system is calcu- 
lated from the resulting density profile n(z) via 



N = n(z)dz. 



(4.12) 



Using the dimensionless coordinate £ and Eq. (|4.11|) 
this can be rewritten as 



N 



Rxmg 



+00 



7(0) 



so that the dimensionless ratio Tq /T is given by 



(4.13) 



T 



2 



j(Rt0 ' 



(4.14) 



This gives a relationship between the global and local 
dimensionless parameters and allows us to present the 
final results in a scalable fashion, rather than in terms of 
the absolute values of N, T, uj z , and g. Here, the desired 
values of Tq /T can be achieved by varying the ratio Ho/T 
of the input parameters (1q and T . 
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FIG. 3: Diagram of different regimes of a uniform ID Bose 
gas in the (7-t) plane. The labels TG, GP, DQ and DC refer 
to the Tonks-Girardeau, Gross-Pitaevskii, decoherent quan- 
tum and decoherent classical regimes, respectively. Although 
all transitions are continuous, for purposes of discussion we 
classify the distinct regimes as follows: 
TG: 7 > 1, t < 7~ 2 ; TG(high-T): 7 " 2 < t < 1; 
GPa: 7 < 1, t < 7" 1 ; GPb: 7 ~ x < t < 7 " 3/2 ; 

" 2 - DC: t > max{l,j- 2 }. 



DQ: 7 " 3/2 < t < j- 



V. DENSITY PROFILES AND PAIR 
CORRELATIONS 

A. Regimes in a uniform gas 

In order to understand the results for the pair correla- 
tions g( 2 >{z,z) of a trapped ID Bose gas, we first recall 
the classification of the regimes of a uniform gas. In Ref. 
[l6|. these were identified using the results for the local 
pair correlation g^ in terms of the interaction parameter 
7 and the reduced temperature r. Here, we give a brief 
summary of these results, except that we rewrite them in 
terms of the parameters 7 and t = T/7 2 , instead of 7 and 
r. This is completely equivalent to the original pair. The 
new parameters 7 and t are more suitable for exploring 
the properties of the trapped g discussed in Sec. 

IV A. 

The diagram representing these different regimes for 
a uniform ID Bose gas in the parameter space (7-^) is 
shown in Fig. [31 The regimes are classified as follows: 

Strong coupling regime: In the strong coupling TG 
regime of 'fermionization', where 7 » 1 and the tem- 
perature T <C Td (t <C 1 or t <C 7~ 2 ), we only have a 
small correction compared to the zero temperature result 

dial: 



7 



4 f 2 



4tt 2 



t < 7 2 , 7 > 1; 



TG. 
(5.1) 



In the case of strongly interacting nondegenerate 
bosons, where 7 3> 1 and the temperature T ^ Td 
(1 <C t <C 7 2 or 7~ 2 <C t <C 1), we have the regime 
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of high-temperature 'fermionization'. Despite the tem- 
perature T » Tj, the loc al p air correlation is strongly 
suppressed (g^ < 1) jHEa: 



.(2) 



~ 2t; j- 2 < t < 1; TG (high-T). (5.2) 



Gross- Pitaevskii regime: In the GP regime, where 7 <C 
1, the chemical potential is fj, = ng and at temperatures 
T <C \i = 2Td"f (t <C 7 or t <C 7 _1 ) the finite temperature 
correction to the zero temperature result is again very 
small El: 



//-• - 1 - -V7 + ^ f V /2 ; 7 « 1, * « 



7 (2) 



GPa. 



(5.3) 

For T > n = 2T d 7 (t » 7 or i » 7- 1 ), the fi- 
nite temperature correction is the leading one. It is 
important to recognize that the upper bound for the 
GP regime extends only up to temperatures of the or- 
der of T ~ -s/^Td (r ~ <fy or t ~ 7~ 3/2 ), and not to 
T ~ Td- Here, the temperature Td is responsible for 
the presence of the quantum degeneracy, while ^fjTd for 
the presence of phase coherence. Thus, for T \x, the 
finite-temperature GP regime lies within the tempera- 
ture interval 2T d 7 < T < ^T d (7 « r « or 
7 _1 <<< 7~ 3 ^ 2 ), and the pair correlation here is given 
by M- 



,(2) 



2 7 ' 



1 « t « 7 _3/2 ; GPb. (5.4) 



Decoherent regime: t 3> max(l, 7~ 3 / 2 ). Due to the 
existence of two characteristic temperatures in the ID 
uniform gas, y^Td and Td, at temperatures T higher 
than ^y^Td one has two sub-regions. For temperatures in 
the interval ^yT d < T < T d (^7 « t « 1 or 7" 3 / 2 < 
i <C 7~ 2 ) the gas is in the decoherent quantum (DQ) 
regime, while for T > Td (t > 1 or t > 7~ 2 ) the gas 
is in the decoherent classical (DC) regime. In both cases 
the local pair correlation is close to g^> ~ 2 |l6j : 



7 (2) 



-,(2) 



2-4/(i 2 7 3 ); 7~ 3/2 «*«7~ 2 ; DQ, (5.5) 



2?r 



1/2 



t > max(7" 2 ,l); DC. (5.6) 



The result in the DC regime remains valid for large 7, 
provided r ^> 7 2 (or £ ^ 7~ 2 ) 16], and we can com- 
bine the required conditions on temperature via t 3> 
max(7~ 2 , 1). 



B. Regimes in a trapped gas 

In a harmonically trapped finite-temperature ID Bose 
gas we again have a strong coupling regime, weak cou- 
pling GP regime, and a decoherent regime. The re- 
sults for the local pair correlation g^(z,z) in the first 



two regimes are easily obtained from Eqs. I|5.1[) - Q5.4|l 
by replacing the interaction parameter 7 by the local z- 
dependent value 7(2) of the trapped sample. However, it 
is convenient to rewrite the results for the local correla- 
tion g( 2 '(0, 0) in the trap center in terms of 7(0) and the 
temperature parameter 9 = T/Tq, where Tq = Nhw z 
is the global temperature of quantum degeneracy of the 
sample as a whole. 

Strong coupling regime. In the strong coupling TG 
regime, where 7(0) 3> 1 and T <C Tq, the density profile 
is given by the Thomas- Fermi result, Eq. (|3.17[) . and 
this allows one to establish the relationship between the 
temperature parameters t and 9, using r(0) ~ tt 2 T/Tq 
valid in this regime. Thus, t ~ tt 2 9/j 2 (0) and Eq. I|5.1|) 
transforms into: 



5 (2) (0,0) 



4tt 2 



37 2 (fj) 



1 



(5.7) 



where 6 < I/71- 2 , and 7(0) > 1. 

For the regime of high-temperature fermionization at 
T ^ Tq, the density profile is given by the thermal Gaus- 
sian, Eq. (g21), so that r(0) ~ 4tt(T/T q ) 2 and hence 
t ~ 47r# 2 /7 2 (0). Therefore, Eq. 1)5. 2jl transforms into: 



5 (2) (0,0H 



8tt 



(5i 



where 1< 6 < 7(0)/v4tt. 

Gross- Pitaevskii regime. In the GP regime (7(0) <C 
1), for temperatures T <C Tq, the density profile 
is given by the Thomas-Fermi inverted parabola, Eq. 
(|3~T31) . We then have r(0) ~ % X /2 1 {Q)T / {ZT Q ) or t ~ 
8%/26»[ 7 (0)]" 3/2 /3, so that Eq. JESJ transforms into 



<? (2) (o,o) = i--Vt(o) + 



16tt 



27 v^Oj 



(5.9) 



where (9 < 3Vt(0)/(8v / 2) (9 < 0.27a/t(0)) and 7(0) < 
1. Similarly, Eq. I|5.4|l transforms into 



<? (2) (0,0) = 1 + 



4V2 
3 



(5.10) 



where 3 v / 7 t0)/(8v / 2) < < 3^/8 (0.27^7(0) < < 
0.27). 

It is important to emphasize that for small 7(0) and 
making the temperature sufficiently low (T/Tq = 9 <C 1) 
the gas is always in the Gross-Pitaevskii regime. 

Decoherent regime. For T 3> Tq, a harmonically 
trapped ID Bose gas is in the decoherent classical regime. 
Here, the density profile is given by the thermal Gaus- 
sian, Eq. 1)4.4(1 . so that r(0) ~ A-k(T/Tq) 2 and hence 
t ~ 47r# 2 /7 2 (0). Accordingly, Eq. I|5.6|l transforms into: 



9 



(2) _ 



= 2- 



7(0) . 
V2 



(5.11) 



where 3> max{l, 7(0)}, neglecting a numerical factor of 
the order of one. Thus, the validity of this result requires 
9 > 1 for 7(0) < 1, and 9 > 7 (0) for 7(0) > 1. 
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FIG. 4: Density profile n(z)/n(0) and the local pair correla- 
tion (z, z) as a function of the distance from the trap cen- 
ter z/Rtf, forT = 0.2T Q (t = 5xl0 4 ) and 7(0) = 1.14xl0~ 3 . 
This corresponds to case lb in the diagram of Fig. HJ below. 



In the case of small 7(0) and at temperatures T <Tq, 
one has a crossover from the classical decoherent (T ^S> 
Tq) regime to the GP (T <C Tq) regime. The proper- 
ties of the gas in this region can be treated as containing 
(locally) features of the decoherent quantum regime of 
the uniform gas. An example illustrating this behaviour 
is given in Fig. 0] where we plot the density profile n(z) 
and the local pair correlation g^(z,z) as a function of 
z/Rtf where Rtf is the Thomas- Fermi radius in the GP 
regime. These are calculated numerically using the solu- 
tion to the Yang- Yang integral equations, with a value of 
7(0) = 1.14 x 10~ 3 at the trap center. 

The temperature in this example is T = 0.2Tq, which 
is intermediate between the decoherent classical and the 
GP regimes. Locally, the tails of the density profile are 
in the decoherent classical regime. On the other hand, 
the central part has features of the decoherent quantum 
regime. The figure for g( 2 >(z,z) shows that fluctuations 
well above the coherent level, with g^(0, 0) ~ 1.5, can 
occur even at temperatures below the transition to a 
quantum gas. However, with further temperature reduc- 
tion well below Tq, the density profile shrinks and one 
always has a coherent GP regime. For smaller values of 
7, the temperature at which the coherent regime emerges 
must in general be calculated numerically. 
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FIG. 5: Local pair correlation g^ 2 ' as a function of the re- 
duced temperature t, for different values of the interaction 
parameter 7. The solid lines are the exact numerical results, 
while the dashed lines correspond to the approximate analytic 
result of Eq. CT . 

C. Variation with temperature 

In Fig. we illustrate different regimes by plotting the 
local pair correlation g^ as a function of the temperature 
parameter t, for different values of 7. 

For sufficiently large 7 the pair correlation approaches 
a universal function of the parameter t: 

v (2) / 2t ' i<cl < (K10\ 

^tJ ] - {2-^, t»i. (5 - 12) 

This is because by increasing 7 one can always reach lo- 
cally the condition t ^S> 7~ 2 . Then, for t <C 1 one has (lo- 
cally) the regime of high-temperature fcrmionization and 
can use Eq. (|5.2I) , whereas for t ^> 1 the sample will be in 
the decoherent classical regime described by Ea. (|5.6() . 

This has an interesting consequence at sufficiently low 
temperatures i <C 1 For 7 > 1, which is always the 
case for far tails of the density distribution, the pair 
correlation remains suppressed below the coherent level 
(g( 2 ) < 1) rather than approaches the value of g^ ~ 2. 
This occurs despite the gas is locally not quantum degen- 
erate at low density. One thus sees that fermionization in 
which the Bose gas develops antibunching with g^ — > 
is an explicitly low-temperature phenomenon, when the 
temperature is scaled relative to the interaction strength. 

However, for 7> 1 the suppression of pair correlations 
is not temperature-independent. Instead, the numeri- 
cal results for increasing 7 converge to a single universal 
function of t. 



D. Spatial variation 

In the case of a trapped gas, the same diagram of Fig. 
|3]also describes the spatial variation of the gas within the 
LDA. The parameter 7 now becomes position dependent, 
j(z), due to the dependence on the density n(z). In this 
diagram, any point on the (7-t) plane can be thought of 
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FIG. 6: Same as in Fig. [3J except with four horizontal lines at 
different temperatures t. This explores different density pro- 
files in the parameter space [7(0)-t], where the points la-lc, 
2a-2c, 3a-3c, and 4a-4c (marked by circles) are representative 
examples corresponding to different values of 7(0) at the trap 
center, at different temperatures t. 



as representing the interaction parameter 7(0) evaluated 
at the trap center, and the dimensionless temperature t. 
This is sufficient to completely characterize the properties 
of the trapped gas. The subsequent local values of j(z) 
of such a gas - as one moves from the trap center towards 
the tails of the density distribution - can be represented 
by a horizontal line drawn in the direction of increasing 
7(2:) at constant t. 

This is shown in Fig. where the four horizontal lines 
correspond to four different temperatures i, while various 
points along each line represent different 'initial' values 
of the interaction parameter 7(0). The interval in the 
left lower corner of the diagram shows the displacement 
(notice the logarithmic scale) for which the local value 
of 7(2;) is increased by a factor of 10. This corresponds 
to a 10 fold decrease in the density n(z). For any given 
distribution with the value of 7(0) in the center, this 
interval helps to immediately determine what fraction of 
the density profile relative to the peak density n(0) is 
contained within a certain regime. 

In the very far tails (z — > 00) of any density distri- 
bution, where n(z) vanishes and 7(2) — + 00, we always 
enter either the DC or the high-temperature TG regime, 
depending on the temperature t. In addition, by con- 
sidering a sample at any fixed temperature t, while the 
peak density is increased [7(0) is decreased], one can al- 
ways reach the situation where the bulk of the density 
distribution is in the GP regime where g^(z,z) ~ 1. 
Physically, this can be achieved by adding more particles 
to the system while maintaining the same global temper- 
ature T, under constant coupling g. From Fig. [S] it is 
clear that the density required may be relatively high, 
with 7(0) = 0.01 being necessary to have a limiting value 
of g^ > 0.9 at t = 10 2 , as an example. For 7(0) > 1, 
there is no coherent GP regime over the entire range of 
temperatures. 




FIG. 7: Density profiles n(z)/n(0) and the local pair corre- 
lation g( 2 \z,z) as a function of z/Rtf for a harmonically 
trapped ID Bose at different temperatures t: t — 5 x 10 4 
(first column) and t — 10 2 (second column). The values of 
the interaction parameter 7(0) in the trap center for each 
of the curves are as follows: la - 7(0) = 1.57 x 10 -4 , lb - 
7(0) = 1.14xl0" 3 , lc- 7 (0) = 0.196, 2a -7(0) = 1.65 xlO -3 , 
2b - 7(0) = 2.31 x 10~ 2 , and 2c - 7(0) = 6.30. The 
dashed lines represent the Thomas-Fermi inverted parabola, 
Eq. 13.131 . The density profiles corresponding to the lower 
density cases lc, and 2c are well approximated by the ther- 
mal (Gaussian) distribution for a classical ideal gas, Eq. 14.41 . 
and are omitted from the graphs for clarity. The respective 
pair correlations for these low density cases are almost con- 
stant along the entire sample and are given by the value of 
g( 2 \z,z) in the tails of the distribution z — > 00. Depending 
on the temperature t, these values can be determined using 
the results of Fig. [SJat 7(2:) 2> 1 [see the curve for 7(2) = 50]. 



To illustrate different examples, we now calculate the 
density profiles n(z)/n(0) and the local pair correla- 
tions g( 2 \z,z) as a function of the distance from the 
trap center z. The distance z is conveniently plotted 
in units of the Thomas-Fermi radius in the GP regime, 
Rtf, given by by Eq. (|3.15(l . The relationship between 
z/Rtf and the dimensionless coordinate £ = z/Rt is: 
z/Rtf = (\/2/4) 1 / 3 t 1 / 2 7(0) 1 / 2 ^. For a gas with a given 
coupling constant <?, the temperature parameter t gives 
the measure of the absolute temperature T. In this sense, 
the examples with t > 1 and t < 1 in Fig. represent 
high- and low-temperature limits, which we analyze sep- 
arately. 



1. High-temperature case 

In Fig. [7| we present examples of calculated density 
profiles n(z) and local pair correlations g^ 2 \z, z), for the 
high-temperature cases of t = 5 x 10 4 and t = 10 2 . The 
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examples shown correspond to the points marked by cir- 
cles la-lc and 2a-2c in the diagram of Fig. HJ1 For each 
temperature t, the sequence of points la, lb and lc (and 
similarly for 2a, 2b and 2c) corresponds to a decreasing 
peak density of the gas [increasing values of 7(0)], while 
the absolute temperature T is kept constant. This can 
be achieved by decreasing the total number of particles 
TV in the sample, at constant T . 

The examples lc and 2c represent a low-density (non- 
degenerate) gas in the decoherent classical (DC) regime. 
The corresponding density profiles are well approximated 
by a thermal Gaussian, Eq. I)4.4|) . and are omitted from 
the graphs for clarity. The respective second-order cor- 
relation functions g^ 2 '(z, z), display large thermal (Gaus- 
sian) density fluctuations with g^(z, z) ~ 2. 

Moving along the horizontal lines in the direction of de- 
creasing 7(0) (starting from the points lc or 2c for each 
temperature t) corresponds to increasing peak densities 
of the gas. As a result one crosses the respective bound- 
aries and enters different regimes of quantum degeneracy 
shown in Fig. [3J Here, the limiting regime as 7(0) — * at 
constant t is the GP regime where the density profiles are 
well approximated by the Thomas-Fermi parabola (|3.13|1 
(see the graphs corresponding to la and 2a), while the 
pair correlation in the bulk of the density profile is close 
to that of the coherent level g^ 2 '(z, z) ~ 1. 

The intermediate values of 7(0) are represented by the 
examples lb and 2b which have density profiles that 
are intermediate between a Gaussian and the inverted 
parabola. The respective pair correlations g( 2 '(z,z) also 
take intermediate values 1 < g^(z, z) < 2. In the exam- 
ple 2b, however, the central part of the density profile is 
in the GP regime, so that the departures form the coher- 
ent level of fluctuations are only seen in the tails of the 
density profiles. 



2. Low-temperature case 

Next, we consider the low-temperature behaviour, in 
which evidence for the Tonks-Girardeau 'fcrmionization' 
can occur. Fig. [S] represents examples of the den- 
sity profiles and pair correlations for a gas with lower 
values of the temperature parameter t: t = 6 x 10~ 2 
and t — 6 x 10~ 4 . The examples shown correspond 
to the points 3a-3c and 4a-4c in the diagram of Fig. 
E| As wee see, for 7(0) > 1 the gas is in the Tonks- 
Girardeau regime. Comparing this with the earlier high- 
temperature examples, we see that for a given density 
with 7(0) > 1, achieving the Tonks-Girardeau regime re- 
quires lower temperatures, t < 1. 

Here again, the low-density [large 7(0)] examples of 3c 
and 4c have density profiles that are well approximated 
by the thermal Gaussian, Eq. ( 14.411 , and are omitted 
from the graphs for clarity. However, the pair correla- 
tions do not display large thermal fluctuations, but rather 
are suppressed below the coherent level, g( 2 \z,z) < 1. 
This reflects the fact that the gas is in the regime of 




FIG. 8: Same as in Fig. except for: ( = 6x 1CT 2 (first 
column) and t — 6 x 10 -4 (second column). The values of the 
interaction parameter 7(0) in the trap center for each of the 
curves are as follows: 3a - 7(0) = 0.323, 3b - 7(0) = 2.16, 3c 
-7(0) = 2.58xl0 2 ,4a- 7 (0) = 5.13, 4b -7(0) = 30.0, and 4c 
- 7(0) = 4.35 x 10 2 . The dashed lines represent the Thomas- 
Fermi inverted parabola in the GP regime, Eq. while 
the dotted lines correspond to the Thomas-Fermi square root 
of parabola in the TG regime, Eq. I3.17H . 



high-temperature "fcrmionization" . The example 4c cor- 
responds to the lowest temperature parameter t, which 
at constant peak density [or constant 7(0)] corresponds 
to the largest interaction strength g, hence the smallest 
value of g^(z, z). 

The example 4b is deep enough in the TG regime, and 
we see that the density profile is close to the Thomas- 
Fermi square root of parabola, Eq. (|3.17l) , while the pair 
correlation is well below the coherent level g^ 2 \z, z) « 1. 
The example 4a is for a smaller value of 7(0) [higher 
peak density], which is closer to the boundary with the 
GP regime. As a result, the shape of the density profile 
departs from the respective TG result and is intermedi- 
ate between the TG and GP parabolas, while the pair 
correlation in the central part of the density distribution 
increases. Finally, the examples 3b and 3a are for inter- 
mediate values of t and 7(0) which are not well described 
by analytical approaches. 

In all these examples the limiting behaviour of the 
pair correlation g^(z,z) in the far tails of the density 
distribution is described by a universal function of t, as 
discussed earlier (see Fig. |SJ. The overall conclusion 
that can be drawn from this analysis is that the local 
pair correlation g^ (z, z) can vary between a broad range 
of values between zero and two and has a rich built-in 
structure. It provides far more sensitive information 
about the regimes of trapped ID Bose gases than the 
respective density profiles. 
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VI. TRAPPED GAS AT CONSTANT N 

Here we investigate the properties of a trapped gas at 
different temperatures T and constant total number of 
particles N. Since the overall picture in terms of the 
density profiles and the behaviour of the local pair corre- 
lation has already been understood in terms of the dia- 
gram of Figs. ElandEl it is now sufficient to only monitor 
the changes in the temperature parameter t and the value 
of the interaction parameter 7(0) under conditions when 
N is kept constant, and then map these changes into the 
[i-7(0)] plane. 

Thus, for a given system with the coupling g, trap 
frequency uj z , and the total number of particles AT, our 
task consists of calculating the density profiles n(z) at 
different temperatures T, with the constraint that the 
total number of particles remains unchanged. Once this 
is done, we identify the respective values of the dimen- 
sionless temperature parameter t and the local value of 
7(0) and plot these on the [i-7(0)] plane of Fig. 

More specifically, instead of performing this analysis 
for absolute values of physical parameters, we first iden- 
tify new dimcnsionless variables for the temperature and 
interaction strength that are more suitable under these 
conditions. The new parameters we introduce are the 
global interaction parameter 7 and the global reduced 
temperature 9: 



Nhuj 7 



= T/Tq, 



1/2 



(6.1) 



(6.2) 



where Tq = Nhuj z . 

The definition of the global interaction parameter 7, 
Eq. ( EH)) ) relies on the identity 



7 2 



r(0) 
7 2 (0) 



= t. 



(6.3) 



Using the definitions of the local parameters t(0) and 
7(0), we see that 7 is the square root of the ratio of two 
energy scales, ma 2 / (2S 2 ) and Tq = Nftui z , as in Eq. (|6.1|) 
(see also Ref. [33|h 

The definitions of the dimensionless temperature and 
interaction parameters 9 and 7 both include the total 
number of particles N. This is more suitable for analyz- 
ing the properties of the gas under conditions of chang- 
ing temperature at constant N. In Fig. Elwe present the 
results of calculation of the density profiles for four dif- 
ferent (fixed) values of the global interaction parameter 7 
while the temperature 9 is changed within a broad range 
of values, typically between 0.1 < 6 < 10. The results 
are summarized by plotting the path of the resulting lo- 
cal values of the interaction parameter 7(0) at the trap 
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FIG. 9: Diagram of the regimes of a trapped ID Bose gas 
as in Fig. except that the curved lines (1) — (4) represent 
the locations of the interaction parameter 7(0) and the re- 
duced temperature t, for four different (fixed) values of the 
global interaction parameter 7 while the global temperature 8 
is changing. This represents four different samples with fixed 
total number of particles N and varying absolute temperature 
T. For each point on a given line, there exists an associated 
density profile with the peak density n(0) corresponding to 
the respective value of 7(0), and the local values n(z) corre- 
sponding to the values of 7(2) in the horizontal direction to 
the right. The values of the global interaction parameter 7 
for each line are: (1) - 7 = 0.002, (2) - 7 = 0.01, (3) - 7 = 1, 
and (4) - 7 = 10. 
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FIG. 10: Variation of the peak density n(0) oc 1/7(0) as a 
function of the temperature 8 — T '/jVftoj z at constant total 
number of particles N. This is the same data as in Fig. |5] 
except plotted in the 1/7(0)-$ plane, where 8 = j 2 t. 



center and the reduced temperature t, in the parameter 
space of Fig. |3| 

For quantitative purposes, we also present the same 
data in the 0-1/7(0) plane, which is shown in Fig. ^| 
While Fig. [5] identifies the local regimes of the gas with 
constant N, Fig. HUI helps to understand the properties 
of the gas in terms of the variation in the global temper- 
ature parameter 9. Note that when 7 is kept constant, 
the variations in the temperature parameters 9 and t are 
essentially equivalent and scale as 9 — j 2 t, according to 
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Eq. 

There are simple approximate relations between the 
global and local interaction parameters 7 and 7(0) at 
high and low temperatures. At high temperatures {6 3> 
1) the relationship is given by: 



7 



7(0) 



[0»1]. 



(6.4) 



At low temperatures (9 <?C 1), and in the limiting GP 
[7(0) < 1] and TG [7(0) > 1] regimes, the relationship 
between 7 and 7(0) becomes independent of 9 and is 
given, respectively, by: 



7 



( 3 



V8\/2 
7- -7(0), 



1/2 



7 (0) 3 / 4 , 7 (0) «1], (6.5) 



[9 < 1, 7(0) > 1] 



(6.6) 



This implies that the GP and TG regimes can equiv- 
alently be defined via 7 or 7(0). The GP regime corre- 
sponds to 7 <C 1 or 7(0) <C 1, while the TG regime will 
correspond to 7 ^> 1 or 7(0) ^> 1. 

From Fig. ^| we see that at high temperatures and 
constant 7, the local interaction parameter 7(0) varies 
according to 7(0) oc 1 / 2 , i.e. linearly in the logarith- 
mic scale of Fig. ED an d in agreement with Eq. ( I6.4fl . 
This means that the response of the peak density n(0) 
to temperature changes at constant N follows the power 
law of n(0) oc T -1 / 2 , which is an expected result for the 
thermal distribution of a classical ideal gas, Eq. (|4.4|) . 

As the temperature is reduced, the response of n(0) 
to the temperature changes becomes modified, and the 
modification is quite different depending on the interac- 
tion strength 7. For weak interactions (7 <C 1), as the 
temperature T is reduced below Tq = Nfku z (9 = 1), the 
peak density n(0) first increases more rapidly than in a 
thermal gas, and then the growth is saturated as the tem- 
perature is reduced further (see curves (1) and (2) in Fig. 
IIOII . At very low temperatures (6 <C 1), the peak density 
n(0) approaches a constant value independent of tem- 
perature. This is a typical behaviour found in a weakly 
interacting gas that undergoes quasi-condensation and 
reaches the GP regime. For intermediate and strong in- 
teractions (7 > 1), on the other hand, the response of 
n(0) to the temperature reduction is different. Instead 
of an initial speed up, the growth of the peak density 
n(0) directly goes to the regime of saturation, once the 
temperature is reduced below the global temperature of 
quantum degeneracy Tq (see curves (3) and (4)). At 
very low temperatures, n(0) again approaches a constant 
value independent of the temperature and the gas ends 
up in the TG regime. 

From the paths of the curves (3) and (4) in Fig. ^\ we 
see that achieving the TG regime from a high tempera- 
ture classical gas by means of reducing the temperature 
T at constant N requires large values of 7 in the first 
place. This can be achieved by having a relatively small 
total number of particles N or a small trap frequency u) z , 
according to Eq. (|6.1|) . 



VII. EXPERIMENTAL CONSIDERATIONS 
A. Average pair correlation 

While the pair correlation g^ 2 \z, z) provides detailed 
information about the local correlation properties of a 
trapped gas, its measurement as a function of z may not 
be an easy task in practice. Here, one usually probes 
the pair correlation (z) 1 ^^ (z) 1 ^ (z) 1 ^ (z)^ within a fi- 
nite volume, e.g. via the measurement of the rates of two- 
body inelastic processes within the entire sample. This 
means that one probes the integrated or averaged corre- 
lation properties of the g MS, S-S has been demonstrated in 
a recent experiment of Ref . . 

We are therefore motivated to study the average pair 
correlation defined via: 



G( 2 ) = I dz(¥(z)¥(z)4>(z)i>(z)) 

dzg^{z,z)n 2 (z). (7.1) 



Here, an interesting question arises of whether this av- 
erage correlation has a simple relationship with the local 
pair correlation at the trap center g( 2 )(0,0). The reason 
to expect this is the fact that g( 2 \z,z) under the inte- 
gral in Eq. i|7.1f) is multiplied by n 2 (z) which vanishes 
rapidly as one approaches the tails of the density pro- 
file. The function g^ 2 \z, z) near the trap center, on the 
other hand, varies slowly and can be approximated by 
the value of g^(0, 0). Therefore, we can approximate 
g^ 2 \z, z) under the integral by a constant g^(0 : 0), thus 
reducing Eq. I|7.1|l to 



G(2) ~ 5 (2) (0,0) / dzn 2 (z) 



(7.2) 



Thus, the average pair correlation can be ex- 

pressed via the local pair correlation </ 2 )(0, 0) using a 
simple relationship, Eq. i|7.2|) . Note that this also re- 
quires an independent evaluation of the integral of the 
squared density, J dzn 2 (z). Introducing a normalized av- 
erage pair correlation, <?( 2 ) , we obtain 



7 (2) 



G( 2 



J dzn 2 (z) 



5 (2) (0,0). 



(7.3) 



In Fig. ^2 we plot the local pair correlation at the trap 
center g^(0, 0) and the normalized average pair correla- 
tion g( 2 ) as a function of the interaction parameter 7(0), 
for four different temperatures t. Each line monitors the 
values of g^(0,0) and g^ as one moves along the lines 
of constant t in Fig. [U Here, the sequence of points 
along the horizontal axis refers to the value of 7(0) of 
the associated density profile n(z), for which we first cal- 
culate the pair correlation g^(z, z) as a function of z 
[which includes the plotted values of g^(0,0)} and then 
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FIG. 11: The local pair correlation at the trap center (0, 0) 
(solid lines) and the normalized average pair correlation g( 2 ~> 
(dashed lines) as a function of the interaction parameter 7(0), 
for four different temperatures t. For each temperature, the 
respective lines monitor the values of (0, 0) and g^ as one 
moves along the lines of constant t in Fig. |^| 

evaluate the integral in Eq. 1|7.1|) to obtain the average 
correlation G( 2 \ and hence g^> . As we see, in the limit 
of small 7(0) the pair correlation approaches the coher- 
ent level of fluctuations with g( 2 \0,0) ~ g( 2 ) ~ 1, while 
at large 7(0) it can take any value between zero and two, 
depending on the temperature t. 

By comparing the full and dashed lines in Fig. 1111 we 
see that the normalized average pair correlation <?( 2 ) can 
indeed be well approximated by the local pair correlation 
in the trap center </ 2 )(0,0). This is an important result 
and may have useful applications in practice. For exam- 
ple, it gives a direct justification of the analysis performed 
in Ref. Q where the results of the measurements of three- 
body recombination rates in a bulk trapped sample have 
been compared with theoretical predictions 0, [nj for a 
uniform gas. 



B. Practical example 

Here, we return to the analysis of Section V, with the 
reference to Fig. and complete it by providing the re- 
sults of calculation of the total number of particles N. 
More specifically, we give the results for the dimension- 
less ratio Nfrui z /T as a function of the local interaction 
parameter 7(0) taken along the horizontal lines of Fig. 

i.e. at four different (fixed) values of the temperature 
parameter t. This is shown in Fig. 1121 For a given trap 
frequency cj z and coupling g, each line corresponds to 
monitoring the variation in the total number of particles 
N as a function of the peak density n(0), at constant 
temperature T. 

Fig. ^]can also be viewed as giving the variation of the 
interaction parameter at the trap center 7(0) as a func- 
tion of Nhu) z /T, which in turn corresponds to monitoring 
the change in the peak density of the gas n(0) as the to- 



FIG. 12: Variation of Nftio z /T as a function of 7(0), for four 
different values of the temperature parameter t. For a given 
coupling g and a fixed T, this monitors the variation of Nhw z 
as one moves along the respective horizontal line in Fig. |S| 
Here, each point along the horizontal line is being referred to 
the value of 7(0) of the associated density profile, for which 
we first calculate the density n(z) as a function of z and then 
evaluate the resulting total number of particles N — J n(z)dz 
to form the dimensionless ratio Nhuj z /T. 



tal number of particles N is varied at constant tempera- 
ture T. Starting from the regime of low particle numbers 
(Nhuj z /T « 1 or high temperatures 6 = T/Nhu z > 1), 
we see that the increase in TV results in a linear increase 
of the peak density, n(0) oc N. This is an expected result 
for the thermal density distribution of a classical ideal 
gas, Eq. (|4.4|) . and corresponds to the linear dependence 
of Nhw z /T on 7(0) as seen in Fig. U3 

As the number of particles is increased further and the 
ratio Nfku z /T goes through the critical region Nfku z /T 
~ 1 (corresponding to temperatures of the order of the 
global temperature of quantum degeneracy, T ~ Tq) the 
growth of the peak density n(0) speeds up, for the lines 
corresponding to t = 5 x 10 4 and t = 10 2 . This speed 
up is most prominent in the first case corresponding to 
very week interactions, and reflects the fact that the gas 
undergoes quasi-condensate formation where the added 
particles mostly go into the trap center rather than into 
the tails of the density distribution. With further in- 
crease of Nhoj z /T, the growth of the peak density n(0) 
is slower than in the classical gas case and the situation 
is now reversed: the added particles mostly go into the 
tails. Here, the gas is deep in the GP regime and the 
density distribution is given by the inverted parabola, 

Eq. prnji . 

For the cases of t — 6 x 10~ 2 and f = 6x 10 -4 , which 
at constant T correspond to stronger inter-particle inter- 
actions, the change in the slop of the respective curves in 
Fig. 1121 around Nftw z /T ~ 1 represents the fact that the 
growth of the peak density n(0) only slows down when 
the number of particles is increased past this critical re- 
gion. Here, the gas goes first through the TG regime 
where the density profile is given by the square root of 
inverted parabola, Eq. I|3.18fl . an d eventually it enters 
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the GP regime (see Fig. with gradual transformation 
of the shape of the density profile to the GP parabola. 

Apart from providing additional information about the 
properties of trapped gases, Fig. ^| can also serve for 
quantitative analysis relevant to practice. As an example, 
we consider a gas of 87 Rb atoms (m = 1.43 x 10 -25 kg, 
a = 5.3 nm) with the aim of identifying a set of physical 
parameters that correspond to the conditions of the point 
3b in Fig. E| Here, 7(0) = 2.16 and t = 6 x 1CT 2 , 
which is an example of a moderately "fermionized" gas. 
The pair correlation g^(z, z) is below the coherent level 
g^(z,z) <C 1, yet the gas is not in the extreme low- 
temperature TG regime and therefore the density profile 
can not be approximated by the square root of parabola. 
In this sense, this example would be easier to realize in 
practice than the extreme TG regime. 

We first consider the trapping potential with lo z /2tt — 
20 Hz and lu ± /2tt = 80 kHz (l z = 2.42 x 10~ 6 m and 
Z_L = 3.83 x 10~ 8 m). With the 87 Rb scattering length 
of a = 5.3 nm, giving g ~ 2hui±a — 5.62 x 10~ 37 J/m, 
and with the value of 7(0) = 2.16 that we are aiming at, 
this set of parameters results in the peak ID density of 
ra(0) ~ 3.35 x 10 6 m ~ x (and hence n 3D ~ 1.82 x 10 20 
m~ 3 ). At this stage, we can identify that the conditions 
of Eq. 12.4fl for achieving the ID regime are satisfied. 

Next, the aimed value of t = 6 x 10~ 2 , together with 
the value of g found above, gives the required absolute 
temperature T ~ 1.22 x 10~ 31 J (in energy units, or T ~ 
8.84 nK). We next refer to the results of Fig. 1 121 and read 
off the value of the dimensionless ratio Nftw z /T ~ 16.7 
(or equivalently 6 = T/Nhui z ~ 0.06), which corresponds 
to 7(0) = 2.16 on the respective t = 6x 10 -2 line. Finally, 
using the values of u> z and T, we find that the required 
total number of particles here is N ~ 153. 

We note that all of the above parameter values are 
close to the conditions realized in recent experiments 0, 
110. 

VIII. SUMMARY 

In summary, we have obtained predictions for the 
correlations and density profiles of a one-dimensional 



trapped Bose gas at finite temperature. This allows pre- 
vious results for the uniform ID Bose gas to be applied 
to the experimentally relevant case of a harmonic trap. 
The calculations use a local density approximation which 
is asymptotically correct in the limit of a large trap with 
a sufficiently slowly varying trap potential. We find that, 
in this limit, there is a similar classification of different 
coherence regimes as in the uniform case. 

Remarkably, the density variation across the trap does 
not cause a dramatic change in the average correlation 
function compared to the value at the trap center. This 
is because the correlations are found to be relatively uni- 
form in the high density region near the trap center, 
which dominates any trap-averaging measurement. This 
is particularly useful for experiments that measure cor- 
relation functions through averaging a nonlinear interac- 
tion over the length of the trap, which is the simplest 
currently available procedure. 

We expect that direct measurements of g^> that can 
test the predictions of this fundamentally important 
many-body theory will become feasible in the near fu- 
ture. 
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